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Abstract. Let / bo a fixed (holomorphic or Maass) modular cusp form. Let 
Xq be a Dirichlet character mod q. We describe a fast algorithm that computes 
the value L(l/2, / X Xq) up to any specified precision. In the case when q is 
smooth or highly composite integer, the time complexity of the algorithm is 
given by 0(1 + |g|5/6+o(l)). 
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1. Introduction 

Let r = SL(2,Z). Let / be a fixed (holomorphic or Maass) cusp form on r\IHI. 
Let Xq be a Dirichlet character on IjqL. In this paper we consider the problem 
of computing central values of L- function corresponding to / twisted by Xq- The 
main result can be summarized as: 

Theorem 1. Let q,M,N he positive integers such that q = MN, where M < N, 
M = M1M2 such that Mi\N and {M2,N) = 1. Let f be a modular (holomorphic 
or Maass) form on T\Ei, s G H and Xq o, Dirichlet character on Z/gZ. Let 7, e be 
any positive reals. Let 

E ^ min{Af^ +iV,(7}. 

Then we can compute L(s, f x Xq) up to an error of 0{q~^) in time 0{E^^'^). The 
constants involved in O are polynomial in (1 + 7)/e. 

This method gives us a positive time saving if q has a factor less than q^/^. The 
maximum saving of size 0(q^/^) can be obtained if q has a suitable factor of size 
« q^/^ . In particular, we can get a saving of size 0{q^/^) for a "smooth" or "highly 
composite" integer q. Note that for these choices of q, the algorithm is considerably 
faster than the 0{q^^°^^^) complexity "approximate functional equation" based 
algorithms. 
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1.1. Model of computation. We will use the real number (infinite precision) 
model of computation that uses real numbers with error free arithmetic having 
cost as unit cost per operation. An operation here means addition, subtraction, 
division, multiplication, evaluation of logarithm (of a complex number z such that 
I arg(z)| < tt) and exponential of a complex number. 

Our algorithm will work if we work with numbers specified by 0(log q) bits. This 
will at most add a power of log q in the time complexity of the algorithm. We refer 
the readers to [101 Chapter 8] and [H] for more details about the real number model 
of computation. 

1.2. Historical background and applications. The problem of "computing" 
values of the zeta function effectively goes as far back as Riemann. Ricmann used 
the Riemann Siegel formula to compute values of the zeta function and verify 
the Riemann hypothesis for first few zeroes. The Riemann Siegel formula writes 
C(l/2 + ir) as a main sum of length 0(r^/^) plus a small easily "computable" error. 
Subsequent improvements for the rapid evaluation of zeta are considered in |13j . 
[T8] . [9], [8], [22], [2], [E] and[l] el al. The current fastest algorithm for evaluating 
C(l/2 + iT) for a single value of T is due to Hiary (time complexity 0(r'*/"+o(i))^ 
see [8]). 

The next natural problem to consider is computing L{\/2 + iT,Xq)^ where Xq 
is a Dirichlet character modulo an integer q. A 0(T^/'^+°(^)(7^/'^+°(^)) algorithm 
for highly composite q, is given by Hiary in |10| . In this algorithm, the rapid 
computation of L(l/2, Xg) is essentially reduced to the problem of fast computations 
of character sums X^feL^o^ Xq{k)^ for any and for M, a small power of q. In case 
of highly composite g, one exploits highly repetitive nature of Xq to get a fast way 
of computing these character sums. The problem of computing L(l/2,Xg) for an 
"almost prime" q seems to be rather difficult. 

In the case of higher rank L-functions, the analogue of the Riemann Siegel for- 
mula is given by the approximate functional equation. A detailed description of 
"the approximate functional equation" based algorithms is given by Rubinstein in 
[15j . In the case of L- function associated to a modular (holomorphic or Maass) 
form, these algorithms have 0(T^+°'^^^) time complexity. The algorithms for rapid 
computation of the GL(1) L- functions unfortunately do not readily generalize to 
the higher rank cases, due to the complicated main sum in the smooth approxi- 
mate functional equation. A geometric way for computing L{f, 1/2 + iT) in time 
(9(2"7/8+o(i)^^ for a modular form / is given in [25] . 

In this paper we consider the higher rank problem corresponding to [10) , i. e com- 
puting i(l/2, / X Xq)- We give a 0(q^/^+°(^)) complexity algorithm for a "smooth" 
or "highly composite" q. Our algorithm in theorem [1] is the first known improve- 
ment of the approximate functional equation based algorithms in in GL(2) x GL(1) 
setting. 

Computing values of L-functions on the critical line has various applications 
in number theory. It can be used to verify the Generalized Riemann Hypothesis 
numerically. It has also been used to connect the distribution of values of L- 
functions on the critical line to the distributions of eigenvalues of unitary random 
matrices via the recent random matrix theory conjectures. 

The problem of computing L-functions is closely related to the problem of finding 
subconvexity bounds for the L-functions. More generally, improving the "square 
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root of analytic conductor bounds" coming from the approximate functional equa- 
tion is of great interest to analytic number theorists. 

In the present paper, we have only considered the GL2 x GLi case. It will be 
of great interest to generalize the method in this paper to higher rank L-functions. 
Integral representations for GL2 x GLi L-functions in the number field setting is 
given by Sarnak in [231 section 11.4]. Our method thus could also generalize for L- 
functions L(l/2, /x x), in the number field setting. A more interesting problem will 
be to generalize our technique in the GL„ x GL„_i setting, where the subconvexity 
bounds are also not known. 

1.3. Outline of the proof. Our algorithm starts with writing L(l/2,/ x Xq) 
essentially as a sum 

9-1 

E/C?/9 + V'?)x,(i), 

3=0 

on the q*'' Hecke orbit of i. The results in this paper are closely related to the 
work of Venkatesh [531 section 6] , where he used the equidistribution of the points 
{j/'J' + */9 ■ ^ J ^ 9" 1} in r\H, to get subconvexity bounds for the L-functions. 
We however use the fact that for a composite q = MiV, we have the following 
decomposition: 

{j/q + : < j < 9 - 1} = uf^-oH(j + kN)/q + t/q:0<k<M- 1}. 

Each arithmetic progression {{j + kN)/q + i/q : < k < M — 1} can be viewed as 
part of the M^^ Hecke orbit of the point Vj = j/N + i/N. 

Let us consider the case {M,N) = 1. We use "well behaved nature" of Xq on 
these arithmetic progressions to convert the problem into the problem of computing 
sums 

^^(^k,cJ{Vj{k)). 
k 

Here, < Cj < M — 1 is such that j = cj mod M and Vj(k) denote the fc*'' point 
in the M*'' Hecke orbit of the point Vj. Here {a/c.;,0 < kj < M - 1} are 
precomputable constants. 

We then "reduce" the points {vj} to points {xj} in a fixed approximate fun- 
damental domain. Let aj be in SL(2,Z) which maps Vj to Xj. Notice that even 
though in r\SL(2,]R), the Hecke orbits of points Vj and xj are the same sets, the 
action of aj permutes them. This means that Vj{k) = aJ^Xj{k) is not necessarily 
equal to Xj{k). However for each j, Vj{k) = Xj(crj{k)), for some permutation aj of 
the Hecke orbit. The problem is now equivalent to computing 

k k ' 

Apriori, there can be M! different permutations of the Hecke orbit. However 
lemma [2T5l shows that the total different number of permutations of the M*^ Hecke 
orbit due to the right action of SL(2,Z) is at most 0{M^). i.e. given M, there 
exists a set of permutations {/3i, /Sa/s}, such that for any Q < j < N — 1, aj = Pk 
for some < fc < - 1. 

Let {ri, rjv/4} be the functions defined on the M*^ Hecke orbit by 
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where < ji < M — 1 and < j'2 < — 1. This inipHes that every j, there exists 
an integer Ij < such that a^-i^i^^ ^, = ri.{k). 

We next use the fact that Hecke orbits of close enough points in the r\SL(2,M) 
remain close. We then sort the points {xj} into sets Ki, K2, ■■ such that the points 
in each set are very close to one other. The number of such required sets is 0{q^'^). 
We choose a fixed representative from each Kj . 

Let us consider Ki for example. Let S Ki be any point. We now need to 
evaluate Sj^ for all jo G Ki. We use a power series expansion around the point 
Xni(k) to compute f{xj„{k)) for each Xj„ g Ki, to get 

d 

^00 ^^Ck^]»J{k,n,,,,Xn^). 
k=a 

Here d ~ 0(1), Ckjo are precomputable constants and J{k, ri, a;„j) can be computed 
for each < I < — 1 and for each k using 0{M) steps. Thus, if we compute the 
coefficients Ckjo and the sums J{k,ri,Xn^) for all Xj„ G ii'i, for all fc < d and for all 

< ^ < M"^ — 1, then we can compute all the sums Sj for all j G Ki "in parallel" 
in further 0(1) time. We can therefore compute Sj for all j in 0((M^ + NY) time. 

Wc thus get a fast way of numerically computing L(l/2, / x Xq)i up to any given 
precision. 

1.4. Outline of the paper. A brief account of the notations used in this paper is 
given in section [21 The algorithm uses a type of "geometric approximate functional 
equation" . It is discussed in detail in section [S] A detailed proof of theorem [T] is 
given in section |4l Lemmas 14.11 and 12.51 deal with well behaved nature of Hecke 
orbits of nearby points in H and well behaved nature of Xq on the "arithmetic 
progressions" respectively. They are proved in section [5l 

1.5. Acknowledgements. I would like to thank Akshay Venkatesh for suggesting 
me this problem. This project would not have been possible without his support. 
Majority of the work in the paper was carried out during my PhD at the Courant 
Institute and Stanford University, and also during my stay EPFL and MPIM, Bonn. 

1 am very thankful to these institutions for allowing me to visit and carry out my 
research 

2. Notation and preliminaries 
Throughout, let F = SL(2,Z). Let / be a holomorphic or Maass cusp form on 

r\H. 

Let s be a fixed point in H. Throughout, let q be a positive integer. Let Xq be a 
Dirichlet character on TLjqL. Let 7,6 be any given positive numbers, independent 
of q. In practice, 7 will be taken to be 0(1) and e will be a small positive number. 

Given a Dirichlet character Xqi The Gauss sum r(xq) is defined by 

9-1 

(2.1) T{Xq)=Y.Xq(hY{klq). 

k=0 

We will denote the set of nonnegative integers by Z+ and the set of nonncgative 
real numbers by M_|- . 
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We will use the symbol <^ as is standard in analytic number theory: namely, 
v4 <C i? means that there exists a positive constant c such that A < cB. These 
constants will always be independent of the choice of T . 

We will use the following special matrices in SL(2,]R) throughout the paper: 

e{x) will be used to denote exp(27rix). 

In this paper, for simplicity let us assume that / is either holomorphic or an even 
Maass form. The algorithms will be analogous for the odd Maass cusp form case. 
For an even Maass cusp form, we will use the following power series expansion : 



(2.3) f{z) = Y,f{n)Wr{nz). 



n>0 



Here Wr{x + iy) = 2^Kir{2'Ky) cos{2ttx) . The explicit Fourier expansion for the 
holomorphic cusp forms is given by 

(2.4) f{z) = fi^X^^)- 

n>0 

For a cusp form of weight k (for the case of Maass forms k is assumed to be 0) , 
the corresponding twisted L- function is defined by: 

Definition 2.1. L{sjx x,) = EZi ■ 

Given a cusp (Maass or holomorphic) form of weight k on r\H, we will define a 
lift / of / to r\SL(2,R) by / : r\SL(2,K) ^ C such that 

(2.) /,(:;;),., „,.r'/(|±^). 

2.1. Real analytic functions on r\SL(2,M). We will use the same notation as 
pS] for real analytic functions on r\SL(2,M). 

Let X be an element of SL(2,M) and let g be a function on r\SL(2,R), a priori 
g{x) does not make sense but throughout we abuse the notation to define 

g{x) = g{Tx). 

i.e. g{x) simply denotes the value of g at the coset corresponding to x. 
Let (j) be the Iwasawa decomposition given by 

^ : (t, y, 6*) e R X M X R ^ n{t)a{y)K{9). 

Recall that (j) restricted to the set R x R x (— 7r,7r] gives a bijection with SL(2,R). 

Definition 2.2. Given rj > 0, let il,, = {—rj^rf) x {-'i],ri) x {—ri,rj) and U,j = 
</)(ii„) C SL(2,R). 

Let us define the following notion of "derivatives" for smooth functions on 
r\SL(2,R): 
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Definition 2.3. Let g he a function on SL(2,R) and x any point in SL(2,M). We 
define (wherever R.H.S. makes sense) 

d d 
^5(2^) ^ lt=o g{xa{t))\ 

d d 

^5(2;) = ^ |t=o g{xK{t)). 

Sometimes, we will also use di to denote 
Given P = {fii, (32, let us define d^^gix) by 

QPi §13-2 

4 

For /3 as above we will define 
and 

|/3| = |/3i| + |/32| + |/33|. 

We now define the notion of real analyticity as follows: 

A function g on r\SL(2,R) will be called real analytic, if given any point x in 
r\SL(2,R), there exists a positive real number r^ such that g has a power series 
expansion given by 



(2.7) g{xn(t)a{y)K{e)) ^ ^ !l||^t/3iyfe0^3 

^=(/3i,/32,/33)eZ3_ 

for every (t, y, 9) S Ur^. 

Let us use the following notation for the power series expansion. 



Definition 2.4. Let y,x G SL(2, M) and t, y, 9 he such that y — xn{t)a{y)K{9) and 
(/3i,/32,/33)=/3eZ^ define 

{y^xf ^t'^'y^'-9^\ 
Hence we can rewrite the Equation (j2.7p as 

/3=(/3i,/32,/33),^ieZ3_ 

Throughout, we will assume that for a cusp form /, all the derivatives of the 
lift / are bounded uniformly on r\SL(2,M) by 1. In general it can be proved that 
given a cusp form /, there exists R such that ||(?^/||oo ^ i?''^', see [24l section 
8.2]. The case when R> 1 can be dealt with analogously. The assumption that all 
derivatives are bounded by 1, allows the proofs to be marginally simpler. 
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2.2. Hecke orbits. Let L be any positive integer and x e SL(2,M), let 

T{L) = {{m,k) : ra\L,{) < k < L/m} 

and 

A{L,m,k) = ^i^ L/m )■ 

The L*^ Hecke orbit is given by left action of the cosets {rA{L, m, k) : (m, k) g 
T{L)}. In particular, the L^^ Hecke orbit of x is given by {rA(L,m, {m,k) £ 
T{L)}, considered as a subset of r\SL(2, M). The Hecke orbits generalize the notion 
of an "arithmetic progression" on SL(2,M). 

It is well known that the right action of any element a in SL(2, Z) permutes the 
cosets {rA{L,m,k) : {m,k) e T{L)}. i.e 

{rA{L, m, k) : (m, k) e T{L)} = {TA{L, m, k)a : (m, k) & T{L)} 

for any a G SL(2,Z). Let CTq : T{L) — > T{L) be be the permutation defined by 

TA{L, m, k)a — TA{L, <7a{m, k)). 

Given any e > 0, using the fact that the number of divisors of M is at most 
O^M"), we get that the cardinality of T{M) is at most 0(A'/^+'^). A priori there 
are M^+'^l possible permutations on T{L). However it is easy to prove the following 
lemma (see section [5]): 

Lemma 2.5. Let ai and a2 are matrices in SL(2,Z) such that ai = 02 modi. 
Then the corresponding permutations and CTa2 0*^6 equal. In other words, for 
every (m, k) G T{L), (m, k) = [m, k). 

Lemma 12.51 implies that the number of possible permutations of the Hecke or- 
bit {TA{L,m,k)\{m,k) £ T{L)} due to the right action of SL(2,Z) are at most 
|SL(2,Z/LZ)| < L^. 

2.3. Specifying / and Xg- We will follow the same assumptions for input of / 
as in |25| . In particular, we will assume that each value of / (or /) or any of its 
derivative can be computed exactly in time 0(1). 

We will also assume that given any integer ri, can be computed in time 

0(1). It can be easily be shown that for q — MN, using and storing a precomputa- 
tion of size 0{M + N), one can compute any Xg("-)i in further O(logq) steps. The 
time complexity 0{M + N) in precomputation does not change the asymptotics of 
the algorithm. Similarly, allowing \og{q) time for each valuation of Xq only adds a 
multiple of log q to the time complexity of the algorithm, which can be absorbed 
into the exponent e. 

In practice, the Gauss sum T{xq) can be computed rather rapidly. A very simple 
0{AP + N) time complexity algorithm can be found in [24[ section 8.6]. As in the 
previous case, this does not change the asymptotic time complexity of the algorithm. 
It is worth mentioning to the reader that the algorithm in |24| section 8.6] is similar 
to the algorithm in the paper and can be helpful in better understanding of the 
underlying idea behind it. 

^ It can be easily shown that given any x and any fixed 7, one can compute d^J{x) up to 
the error 0[q ') in 0(q°^^') time. Here the constant involved in O is a polynomial in and 
7. In this algorithm we only compute values of d^f{x) for \/3\ <C 1. Allowing 0(g°'^') time for 
each valuation of / docs not change the time complexity of the algorithm. See 1241 chapter 7] for 
explicit details about it 
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3. A GEOMETRIC APPROXIMATE FUNCTIONAL EQUATION FOR L{s, f X Xq) 

Our algorithm will start with proving a "geometric approximate functional equa- 
tion" for L{sJ X Xq), given by ([XT^ and ^A^. The right hand side of ([51^ 
and (|3.13p consists of sums of q integrals. Each of these integrals is an integral of a 
'nice' function on a geodesic of (hyperbolic) length O(logq). Therefore, using [25l 
proposition 8.1], we can write each of these integrals (up to an error of 0{q^^)), as 
a sum of size 0{q'^) terms. Adding all these sums together, the right hand sides of 
p.l2p and p.l3p can be written (up to an error of 0{q~'')) a sum of size 0((7^+^). 
The constants involved in O are polynomial in (1 + "f)/e and are independent of q. 

Holomorphic case: 

q-l 9—1 CXD 

X! Xq{k)f{k/q + iy)^Yl Xq{k) ^ f{n)e{nk/q)e{iny) 

k=Q fe=0 n=l 

oo q — 1 

= Xqik)e{nk/q) 

n=l k=0 
oo 

(3.1) = T(xq) ^ /(n)x9(n)e(m?/). 

n=l 

After taking the Mellin transform of p.l|) , we get 

9-1 j-oc 

(3.2) ^X,(fc) / f{k/q + iy)y^+('^-^y^dy 
Jo 



k=0 



„s+(fc-3)/2 



dy, 



poo 

^(X?) y] /HX9("-) / e.{iny)y' 
n=i 

pOO 

r{Xq)y,f{n)Xq{n) / exp(-27rny)y^+('^-3)/2dy; 
„=i -^0 

^^-^g^|LL(.,/xx,)r(. + (fc-l)/2). 



Even Maass form case: 
9-1 

Y.^q{k)f{k/q + iy) 



k=0 



9 — 1 00 

X! f {n)\/ny cos{2TTnk / q)Kir{2TTny) 

k=0 11=1 

2 f /(n)VH^X.(2.ny) g ^^^^fJl^^^t^ 

n=l k=0 



T 



iXq)+TiXq) 



Y fin)Xq{'n)^/riyK,r{27Tny) 



(3.3) = r(x,)(l + X9(-l)) Y f(n)xq{n)^K,r{2nny). 



n=l 



Taking Mellin transform of p.3p , we get 
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9-1 roo 
(3.4) 5]xg(fc) / I{k/q + iv)v^-^'^dv 
k=o -^0 



= r(x,)(l + X,(-l))E/W^^9(") / K,r{27Tny)y'-^dy; 

n=l 

= -(X.) (2,). ^(^. / X x,)r [-^ ) r 

r(xg) is a complex number with absolute value g'. Therefore we can use p.2p 
to compute L(s, / x x^) in the holomorphic case. 

For an even Maass form / however, if Xg(~l) = then the right hand side 
of p.4p is zero. Therefore it needs slightly different treatment. Recall the Fourier 
expansion for / given by 

f{z) = ^ f{n)2y^K,r{2TTny)cos{2Trnx). 

n>0 

This implies that 

dxf{z) = —2n nf{n)2y/yKir(2TTny)sm(2'!Tnx). 

Tl>0 

We use a similar method as before to get: 

9-1 

^Xqik)dxf{k/q + iy) 

q~l oo 

= —2n Xq{k) ny/riyf{n) sin(2Tmk/q)Kir{2TTny) 

k=0 11=1 

e{nk I q) — e{—nk I q) 



-2^ f{n)n^'^^K,r{2^ny) ^ Xgl^)- 



fc=0 



2i 



n=l 

oo 

(3.5) = z^t(x,)(1 - X,(-l)) E f{nW\,{n)^K,ri2T,ny). 

n=l 

Taking McUin transform of (|3.5I) . wc get 

9-1 poo 

(3.6) ^xg(fc) / 9./(fc/'? + *y)y'''/'dy 

fc=0 -^0 



K„{2TTny)y^dy 



»^r(x,)(l - X,(-l)) E /HX, W"'/' / 
„=i -^0 

--(x9)^^^yiT^i(-, / X X9)r — ^ — j r ^ — 
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Notice that p.6p is analogous to p.4p . The algorithm to compute 
X]fc=o Xg(fc) Jo° 9xf{k/q + iy)y'^~^^^dy is completely analogous to the algorithm 
to compute J2k=o Xgi^^) /o°° f{k/q + iy)y''~^^^dy. Hence throughout the rest of the 
paper, we will assume that x<?(l) = Xq{~^) = 1- Using the automorphy of /, we 
get the following lemma (analogous to [ISl Lemma 3.1]). 

Note that a similar treatment will give us the corresponding "geometric approx- 
imate functional equations" for odd Maass forms. 

Lemma 3.1. Given any cusp form f (of weight k) on r\SL(2, M), positive integers 
q such that n < q, s G H, a positive real 7 and for any c > 2, 

(3.7) / f{n/q + iy)fdy^ / f{nlq + iy)fdy + 0{q-'). 



The constant involved in O is independent of q and c. 

Proof. We begin by noting that it is enough to prove the lemma when n is coprime 
to q. Let us use the exponential decay of / at ioo, to get for y > 1, 

(3.8) |/(n/g + zy)| <exp(-7ry). 

This implies that we can choose a constant ci, independent of n, q such that for 
every y > ci log q, we have 

(3.9) |/(n/g + zy)|«q-l^|-^-V'. 

Let n',n" such that < n' < q and nn' — qn" = 1. The the action of g = 

on H maps n/q to infinity. Using the automorphy of / with respect 



-q n 
to the action of 5, we get 

(3.10) f{n/q + ^y) = {-q{n/q + ^y) + ("^ +^^) ' ) 

—q(n/q + ly) + n 

= i-qw) / ■■ — 

-qyi 

^{-qiy)-'f{-n'/q+^). 

Q y 

We use the exponential decay of / at infinity to get a constant C2 , independent 
of n, g such that for every y < ■, we have 

(3.11) |/(n/g + zy)|«9-^-l^l-^ 
The equations p.9p and (|3.1ip give us the result. 

□ 



Lemma l3. II implies that for any c > 2, 
(3-12) ^^,yiLy. L{s, f X XqMs + (fc - l)/2) 

= Ex,(j) r fij/q + iy)y'+^'-'^^^'dy + 0{q-'^). 

and similarly we get that given any (7,7 > 0, a Maass cusp form /, and c > 2, 
we have 
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(3.13) rU.)'-^^Hs, f X x,)r {^) r (^) 

= Ex5(fc) r /(fc/g + z2;)r'-My + 0(g-^). 

The equations p.l2p and p.l3p denote "geometric approximate functional equa- 
tions" to compute L{s, fxxq) in the holomorphic and Maass form case respectively. 
Notice that the integrals on the right hand side of p.l2p and p.l3p are over hyper- 
bolic curves of length <C logg, therefore they can be computed in 0(g°(^)) time. In 
the following theorem, we discretize the integrals in these equations to convert the 
problem of computing the sum on the right hand side of p.l2p / p.l3p to the problem 
of computing the sum S ~ S|=o Xg(j)f^2/("(j/'i')'^(^)) f^'" ^i^Y ^ ^ [—clogq, clogq]. 
Here 82 denotes jr-- 

Lemma 3.2. Let f he a modular (holomorphic or Maass ) cusp form on T\M., and 
Xq be a Dirichlet character modulo q, s be any complex number. Let 7, e be any 
positive reals, then there exists a positive integer N' = 0((1 -f 7)/e) such that if for 
any \t\ ^ logg and for any < I < N' , we can compute X]j=o Xq{j)d\f{n{j / q)a{t)) 
up to a maximum error 0{q~^) in time D{q), then we can compute T{xq)L[s, f^Xg) 
using 0{D{q)q'^) operations. The constants in O are polynomial in (1 +j)/e. 

Proof. Let us start with the "geometric approximate functional equations" 
p.l2p / p.l3p for holomorphic/Maass case respectively. We use the lift / defined 
in ([23]) to get /(n(x)a(log y)) = y'^^^fix + iy) to rewite (f3J2|) as: 



^^^^^^,L{s,f.Xq)ns + {k-l)/2) 

9-1 /.q" 



r fi^i3/q)ci{\ogyW-'"^dy + 0{q- 



i=0 

Notice that the above equation will be valid for any c > 2. Substitute logy = < in 
the above equation to get that 

^^^p^i,y, L{s,f^xq)n^ + {k-i)m 

= EX,(J) r'°'' f{n{j/q)a{t))e'(^-''^Ut + 0{q-^). 
j=0 -'-clogij 

Therefore we have 

l^'^'J J -clog q 

Here C = = sup_,i„,„<,<,i„_ | exp(t(s - 1/2))| 



-c log g<t<c log g I 

1 



9 At) = 7^/("(j79)a(i))exp(t(s - 1/2)). 
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Notice that ^ \t=to {fin{j/q)a{t))) = ^ f {n{j / q)a{to)) . Here ^ denote the 
derivative of / in the "geodesic direction" , defined in section [2] 

We have assumed that / has bounded derivatives (ref. section [2]). Therefore 
using Leibnitz rule, for a fixed s S C and each t in [— clogg, logg], 

|^5,(to)«/(M- 1/21 + 1)". 

The constant involved is independent of q. Hence for each t in [— clog g, clog g], 
gj is real analytic with radius of convergence at least l/(|s — 1/2| + 1). Therefore, 
choosing a grid of 0{q'^) cquispaccd points and using power series expansion at the 
nearest grid point on the left to compute g}^ at any given point, we get that given 
any 7, e > 0, 

f{n{j/q)a{t))cMt{s - l/2))dt 

'-clogg 

c9'(logg)-l N' q-' ; 

(3.15) =C' J2 Ey diig,{xq-^))-dt + 0{q~'~-') 

x— — cq'^ log <? /— ^ 

Here N' = 0((1 + 7)/e). Notice that the above equation is true for any c > 2. 
Hence , given e, we can choose c such that {cq'^ log q} = 0. Multiplying p.l5p by 
Xq(fc) and summing over fc, we get 

cq^ log q—1 N' q~l 

(3.16) L{sJxxq)^C" £ E^'E^'?(-?)^2fe(-^9-^)) + 0(g-^) 

x — — cq^ log g /— J — 

Here d, - '-dt - and C" - Ci^-r^"-'"' 

° V. (i+iy. 'I'liu - r(s+(fc-i)/2)T(x,)- 

(j3.16p , along with the definition of gj (t) implies that the problem of computing 
L{sJ X x) is equivalent to computing T{xq) and Y.'jZl Xq{j)dif{n{j /q)a{t)) for 
any t € [— clog q, clog q] and for < Z < A^' faster than 0{q). Here 82 denotes 
Exactly same method will work for the even Maass forms, when Xg(— 1) = — 1- 

Note that for an even Maass form if Xq(— 1) = —1, then we need to compute 

the sum JZjZl XgU) Iq-c dxf{j /q + iy)y''^^^'^dy. Notice that f{n{x)a{\ogy)) = 
f{x + iy). Therefore, for any xq, an easy computation gives 

d ~ 

dxf{xo + iy) = j/"^^/(7i(.To)a(logy)). 

Using this equation, we can follow exactly the same method before, to get lemma 
when Xg(~l) = ^1- D 

In the following sections, we prove that D{q) = (A/^^"^ + A^)^+'^. We use this 
result to finish the proof of theorem [1] 



4. Proof of theorem [T] 

Let r ~ SL(2,Z) and / be a (holomorphic or Maass) cusp form of weight k on 
r\EI. Let q = MN where M < N, M = M1M2, where Mi = gcd(M,iV) and 
{M2,N) = 1. Let Xq be a Dirichlet character on "LjqL. Let 7,6 be any given 
positive numbers. The main goal of this section is to prove theorem [T] 
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Lemma l3. 21 implies that the algorithm to compute L(s, / x Xq) faster than 0{q) 
is equivalent to computing the Gauss sum T{xq) and 

q-l 

(4.1) Y.^<i{k)dif{n{k/q)a{t)) 

k=0 

for any t E [— clog q, clog q] and for < I < N' , faster than 0{q). Here 82 denotes 
N' ~ 0((1 +7)/e) and c > 2 is a suitable computable constant. Hence, in this 
section we will consider the problem of computing the sum in (|4.ip for any given 
t eJ-clogg,clogg]. 


Computing T(xq) can be done very rapidly. It is easy to see that given any 
positive 7, we can compute T{xq) up to error 0{q~'''), using at most 0{M'^ + N) 
operations. 

Wc will deal with computing (|4.ip in the special cases: q — AIN where (M, N) ~ 
1 and the case when M\N separately in the subsections 14.11 and 14.41 respectivelv. 
We use these results to complete the proof of theorem [T] in section 14.51 

4.1. Case q = MN where {M,N) = 1. Let g be a real analytic function with 
bounded derivatives on r\SL(2,]R). Let q = MN where {M,N) = 1. Let t be a 
real number in [— clog g, clog (7], ti = t + logq. Let 

xo = a{ti). 

In this section, we will consider the problem of computing the sum 

9-1 

5 = E^9(%KV'?)a(i)) 

fc=0 

9-1 

= Y.Xqik)g{A{q,l,k)aih)) 

k=0 

N-1 M-1 

(4.2) = E E X<i{j + Nk)g{A{q,l,j + kN)xo). 

j=0 k=0 

Computing the sum of this type is equivalent to computing the sum (|4.ip . Recall 
that A{q, 1, fc) is defined in section O Use, A{q, 1, j + kN) = A{M, 1, k)A{N, 1, j) 
to rewrite (|4.2p as 

(4.3) S^Y.^r 

j=o 

Here. Sj is defined as 



Note that for a Maass form /, for g = f, t = — logg and I 
form 

9-1 
fc=U 

In this section, we essentially give an algorithm to compute 
0{q). 



= 0, 114.111 takes a more familiar 



Xq{k)f{k/q + i/q), faster than 
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M-1 

(4.4) S, = + ''N)9iMM, 1, k)v,). 

k=0 

Here 

Vj = ^(7V,l,j>o. 

We now use Xg = XmXn to get 

M-l 

'5,= ^ X,(j + fciV)5(^(M,l,fc)«,) 

fc=0 

A/-1 

(4.5) = xnU) Xa/(j + kN)giA{M, 1, fc)«,)- 

fe=0 

Let us "reduce" the points to an approximate fundamental domain 

using the algorithm in |24[ Chapter 7] . Hence we have matrices {7^ , xj } such that 
Vj = "fjXj, 7j € SL(2,Z) and Xj lies in the "approximate fundamental domain". 

Xm is a character on Z/AfZ. Given j, let Cj be defined by j = Cj mod M and 
< Cj < A/. Hence rewrite (|4.5p as 

A//m-l 

(4.6) 5^=XAr(j)5I I] /ic,(m,%(A(Af,m,fc)7,x,). 

m|A/ k=0 

Here for {0 < ^ < A/ - 1}, hi : T{M) is defined as: 

(4.7) hi{{m,k))=5i{m)xM{l + kN). 

Where, 5i is the characteristic function of {!}. For example, for a prime M the 
functions hj are given explicitly by: hj{\, k) ~ xmU + kN) and hj{M, 0) = 0. 

The number of points in T{M) are at most Af^+' ( using the fact that the 
number of divisors of A/ are at most 0(Af^)). The right action by a € F permutes 
the T{M) in F\SL(2,M). In other words, given any element a of F, there exists a 
permutation Oa on T{M) such that for each {m,k) G T{AI), we have an a' G F 
such that 

A{M,m,k)a = a'A{M, aa{m,k)). 

Using lemma [^751 we get that if a = 6 mod Af, then aa = crt- This implies that 
the total number of permutations of T{M) due to the right action of F is at most 
M^. Using this, we rewrite (|4.6p as 

(4.8) S.^xnU) Y K^{a^^{m,k))g{A(M,m,k)xj). 

{m,k) 6 T{M) 

Lemma [^31 implies that the number of distinct functions hc^cr^. is at most M^. Let 
us enumerate them as ri, ...,ri (say) , where L < Af^. 

Notice that the Hecke orbits "preserve" the distance between the points. In other 
words, let M be any positive integer , and x,y £ SL(2, K.) such that x G yV for some 
open neighbourhood V of the identity. Then A{M, m, k)x lies in A{M, m, k)yV, for 
all {m,k) £ T[M). Therefore, given any positive e, we "sort" the points {xj} into 
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sets Ki,...,Kq such that x,y £ Kj imphes that x^^y G Uq-t. Here Q = O((j'^')0 
We choose a representative Xn^ from each Kj . 

Let us focus on Ki. Let us use a power series expansion around the point 
A{M,m,k)xni to compute g(A{M,m,k)x) for all x G i^i. We summarize the 
result in the following lemma: 

Lemma 4.1. Given any positive reals 6,7, positive integer M, x and y € Ki for 
some integer i. Let r be any complex valued function defined on the the set T{M). 
Then there are explicitly computable constants cp^x,y (^nd d such that, 

d 

(4.9) J(0, r,y)=Y, Cf,,,,yJ{P, r, x) + ©(M^g"^). 

|/3|=0 

i/ere d = 0((1 + 7)/e) and J{p,r,x) is defined by 

(4.10) J(/3,r,x)= r{{m,k))d^g{A{M,m,k)x). 

The constant in O is independent of AI, q and is a polynomial in (1 + 7)/e. 

We will prove the lemma HTT] in section [5] Notice that for fixed /3,r and a fixed 
X, we can compute J(/3, r, x) in time 0(M^+'^). This gives us an exact idea of the 
algorithm. The exact algorithm is given by: 

4.2. Explicit algorithm. 

(1) Compute and store the functions ri, ...,rL. 

(2) Reduction of points vi, ...,vn into points xi, ...,xn and sorting of the points 
into sets Ki, ...,Kq and choose a representative x„. from each Ki. 

(3) For each < ^ < A^, compute < q < M — 1 be such that q = I mod M. 
Find ji such that /ic, cr-yi = rj^ . 

(4) i ^ 1. 

(5) for all xi in Ki, and for all |/3| < d, compute and store C/j^x„.,xi- 

(6) for all |/3| < d, and for all 1 < t < L, compute and store J(/3, rt, 

(7) for each xi in Ki, compute 

d 

Si^Xn(I) Y ^l3,Xr,.,x,J{(3,rj^,Xn,)- 

\I3\=0 

(8) if i = Q compute S = Si + ... + Sn else i ^ i + 1 and go to step [S] 

4.3. Time complexity. Let us compute the time complexity in the algorithm. 
The O constants involved in this proof are polynomial in (1 + 7)777. 

Computing and storing for all i takes at most (roughly) 0{LM^^'^) time. But 
recall that L < M^. Hence step[I]takes 0(Af^+") time. 

It is easy to see that "reducing" each Vi to an approximate fundamental domain, 
can be done in 0(log q) time. There are many standard reduction algorithms avail- 
able. We refer the readers to [55] for a form of the reduction algorithm. Cutting 
the approximate fundamental domain into Q = 0{q'^'^) "cubes" and sorting the 
points Xj in the sets Ki, Kq, and choosing a representative Xm from each Ki 

■^The approximate fundamental domain has form {n{t)a(y)k[6) :0<Ct<Sl,l<C3/<C 
logg,— 7r < 6 < 7r}. We divide it into Q "cubes" having sides 0(g~"^) each. Every Xj lies in 
one of these Q cubes. For each j < Q, Kj consists of points Xi is in the j^^ cube. 
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takes 0(iVi+' + q^^) time . Step [3] also takes 0{N) time. The whole "reducing" 
and sorting process has also been discussed in detail in [24l chapter 7 ]. 

Notice that for each xi G Si, we need to compute Ci3^x„.,xi for <C 1. In 
sectionlSl it will be shown that each of these values can be computed in 0(1) time. 
Therefore, step[5]takes 0(1) steps. Hence for all I, step[5]takes 0{N) steps. Notice 
that for a fixed i, and fixed t, steplHltakcs 0{M^'^'^) time. 

Hence for fixed i, and for all t, step |6] takes 0{LM^~^'^) time. Recall that L < 
. Therefore total time spent in computing J(/3,rt,a;„J for a fixed i and all t is 
0{M^+^). Therefore, for all i, for aU j, and for all i, the step|6]takes 0{QM^+^) w 
0(M5+4^) ^ij^g^ RecaW here that Q « M^^ and d = 0(1). Steps [7] and |8] take 0{N) 
steps. 

Therefore, the total time spent is 0(M^+'^" + N). 

4.4. Case q = MN when M\N . Let q = MN where M\N . Unless redefined in 
this section, we borrow the notations from previous section. 
We proceed as in previous section and rewrite (|4.2p as 

(4.11) ^=E^^- 

Here, Sj is defined as 

M-l 

(4.12) S, = + kN)g{A{M, 1, k)v,). 

k=0 

Vj as in the previous section. For {j,N) 7^ 1, Sj =0. Therefore, rewrite (|4.1ip as 

(4.13) 5, = J2* S,. 

j mod 

Let F{k) = X(j(l + kN). Since F(fci + fc2) = F{ki)F{k2), there exists a constant 
b such that F(fc) = Xq(l + fcA^) = e{bk/M). Using this we have that for any j, 
coprime to N, 

M-l 

(4.14) S, = x,U) E 4bj-^k/M)g{A{M, 1, k)v,). 

k=0 

Notice that denotes the multiplicative inverse of j mod m. Let Cj be defined 
by Cj = j mod M, where < Cj < q — I. We rewrite (|4.14p as 

M/m-l 

(4.15) Sj^xgif) E K^{(m,k))g{A{M,m,k)vj). 

(m,k)eT(M) 

Here, for < ? < M -l,hi: T(M) 

(4.16) hi{{m,k)) ^ Si{m)e{bl-^k/M). 

The function hi only depends on the residue of I mod M. Therefore, there are 
M distinct functions hc^ ■ We can proceed exactly similarly as in last section to get 
the required algorithm. 
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4.5. Proof of theorem [H Let q = MN. Let Af = M1M2 such that {M2,N) 1 
and A/i|iV. We proceed as in previous sections and rewrite ()4.2p as 

(4.17) S^Y.^r 
Here, S'j is defined as 

A/-1 

(4.18) S, = X,{j + kN)g{A{M, 1, fc)«,). 

Vj as in the previous sections. Notice that Xg — XM2XM1N and that there exists b 
such that XMi7v(1 + kN) ~ e{kb/Mi). Moreover, Sj = 0, if j is not coprime to N. 
Therefore for {j,N) ~ 1, rewrite (|4.18p as 

A/1-1 A/2-1 

= E E XgU + (fc + lAh)N)g{AiM, l,k + lAh)v,) 

k=o 1=0 
A/1-1 A/2-1 

= £ £ XqU + kN + ahN)g{A{M,l,k + lAh)v,) 
k=o 1=0 

A/1-1 A/2-1 
= £ XM,NU + kN) J2 XAhU + kN + lAhN)g{A{M,l,k + lAh)v,) 

k=0 1=0 

(4.19) 

A/1-1 L'-lj, A/2-1 

= XA/iJv(j) 51 ^(^f^) E XA/2(i + A:Af + ^MiiV)g(A(Af,l,fc + ?A/i)z;j)- 

Notice that for fixed k and the quantity e{bj^^k/Mi) is uniquely determined for 
j mod Afi and XM2U + kN + IMiN) is uniquely determined for j mod M2. Hence 
if ji = j2 mod Af then for all < fc < A/i - 1 and < Z < Af2 - 1 wc have that, 

<-\i-)XA'h Ui +kN + lAhN) = e(:^\— )xA/2 (j2 + kN + lAhN). 
Ml All 

Let Cj be defined by Cj = j mod Af, where < Cj < AI — 1. We can rewrite (|4.19p 

as 

(4.20) Sj=xm,n{j) K^{{m,k))g{A{AI,m,k)v,). 

Here, for {l,N) = 1, /i; : T(Af) ^ is defined by: 

(4.21) hi{{m,k)) ^ 5i{m)e{br^k/Adi)xMAl + koN + loAIiN)- 
where, 

= fco + IqAIi, < fco < A/i - 1, < ?o < A/2 - 1. 

The inverse in (j4.2ip is modA/i. Notice again that there are only at most AI 
distinct functions he on T{M). Therefore, the algorithm in l4.1l can be easily used 
to get an algorithm in the general case. 
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5. Proof of lemmas [47T] and [2751 
In this section we will give a proof of lemmas 14.11 and 12.51 

5.1. proof of lemma [4. li Let x-, y G Ki for some integer i and let r be a function 
on the set T{M). This implies that x^^y G Uq-c This implies that 

{A{A4,m,k)xy^A{M,m,k)y = x~'^y G Ug^e 

for aU {m,k) G T(Af). Let 

x~^y ^ n{to)a{yo)K{ea). 

Here \6o\ < tt. As 17 is a real analytic function on r\SL(2, R), we can use the power 
series expansion for g at points A{M, m, k)x to compute g{A(M, m, k)y) to get 

(5.1) ,(A(M,„,fc),)= ^^^((-^(^;'-'^)-) ,g^,g-<3. 

As n{to)a{yo)K{9o) G Uq-e, we get that there exists a constant d ~ 0((1 + j)/e) 
such that 

(5.2) yiAiM,m,k)y) . fe).) ^,,^,.^,3 ^ o(,-^). 

l/3|<d ^' 

Substituting in (|4.8p . we get 

d 

(5.3) J(0, r, y) = ^ cp,^,yJ{l3, r, x) + 0{2M^q-^). 

1,91=0 

Here cp^x,y = *o'2/o'^o'7/3!- Notice that each C/3,;,n can be computed in unit time. 

5.2. proof of lemma l2.5i Let (m, fc) be a any element of T{L). Let aa^{m,k) = 
(mi, ki). This implies that TA{L^ m, fc)ai = rA(L, mi, fci). This implies that there 
exist a matrix c G SL(2, Z) such that 

(5.4) cA(L,m,k)ai — A{L,mi,ki). 

The lemma is equivalent is proving that there exists a matrix d in SL(2,Z), such 
that 

(5.5) dA{L,m,k)a2 ~ A{L,mi,ki) 
( this would imply that ("*; k) — (Tai {'m, k)). 
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We use (|5^ to get 

cA(L, TO, k)a2 
~ cA(L, TO, k){ai + Lb) 
= cA(L, TO, k)ai + cA{L, to, k)Lb 
= A{L, nil, ki) + cA{L, to, k)Lb 

= A{L, TOi, fci) + cA{L, TO, k)LbA{L, mi, ki)~^A{L, mi, ki) 
= (/ + cA{L, m, k)LbA{L, toi, toi, fci) 

TO k \ 7 f L/mi —ki V_i 



I^M implies that d = c(/ + c ^ , 6 '„ ^ / will be in 

\^ L/m J y O L/TOi J 

SL(2,Z) and will satisfy condition in equation (j5.5l) . {m,k) was any arbitrary 

element of T{L). This implies that = Oa^. 
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